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I .  INTRODUCTION 


This  report  deals  with  the  problem  of  optimizing  the  packaging  for 
artillery  shells  to  provide  protection  against  influence  initiation  in 
the  metal  jetting  environment  produced  by  shell  casing  impacts.  The 
initiation  of  a  shell  stored  or  transported  with  other  shells  produces 
high  speed  casing  fragments  that  on  impact  with  nearby  shells  can  cause 
shocks,  penetration,  and/or  fracture  resulting  in  the  initiation  of  the 
explosive  fill  in  the  neighboring  shells.  In  addition,  the  detonation  of 
two  adjoining  shells  can,  as  their  casings  impact,  produce  significant 
inetal  jetting,  In  the  case  of  15S  mm  artillery  shells  this  metal  jet  can 
cut  through  as  much  as  15  cm  of  steel.  A3  such,  this  latter  mechanism 
represents  the  most  severe  fratricide  environment  constraining  shield 
design.  Neighboring  shells  can  be  protected  from  this  jet  of  material  by 
inserting  layers  of  material  to  intercept  and  disrupt  the  metal  jet.  The 
insertion  of  material,  however,  greatly  increases  the  total  amount  of  material 
required  for  the  safe  packaging  of  munitions.  It  becomes  important,  there¬ 
fore,  to  develop  a  model  of  the  shield  interaction  with  such  shell  casing 
jets  to  facilitate  the  development  of  optimum  packaging  designs  capable  of 
inhibiting  influence  initiation  of  shells  in  this  severe  environment. 

The  design  of  an  optimum  packaging  involves  variation  of  numerous  para¬ 
meters  such  as  number,  thickness,  composition,  and  placement  of  the  shield 
materials.  Extensive  experimental  work  has  been  carried  out  to  determine 
good  potential  designs1 2 3”4.  Still,  determination  of  an  optimum  design  requires 
the  development  of  a  mathematical  model  to  minimize  the  total  number  of 
experimental  tests  necessary  for  the  development  of  the  influence  suppression 
configuration.  The  purpose  of  the  present  report  is  to  provide  a  model  of 
the  antifratricide  shield  panels  subjected  to  the  most  severe  environment, 
jets  produced  by  simultaneously  detonating  neighboring  rounds.  A  subsequent 
report  will  deal  with  the  problem  of  shock  attenuation  by  panels  separating 
neighboring  rounds. 


1  Howe ,  P.  M, ,  " The  Phenomenology  of  Interround  Communication  and  Techniques 
For  Prevention,"  ARBRL-TR-02048  (1978)  (ADHA054373) . 

2)Iowe,  P,M, ,  Collie,  D. ,  "Effectiveness  of  Plastic  Shields  in  Pi'evention  of 

Propagation  of  Reaction  Between  Compartmentalised  Warheads,  "  ARBRl-MR-0282? 
(1978)  (ADHB027466L). 

3Howe,  P.M.,  "The  Response  of  Munitions  to  Impact, "  ARBRL-TR-02189  (1979) 
(ADHB040230)  . 

4 Gibbons ,  Jr,,  Gould ,  " Multiple  Round  Fragmentation  Hazards  and  Shielding," 
ARBRL-TR-02329  (1981)  .  (AD  U0S8793L) 


II.  DESCRIPTION  OF  THE  ARTILLERY  SHELL  METAL  JET 

The  mechanism  causing  influence  initiation  of  concern  in  the  present 
report  is  the  jetting  of  metal  due  to  the  impact  of  shell  casing  material 
from  two  neighboring  shells  initiated  at  approximately  the  same  time.  An 
example  of  the  formation  of  a  jet  resulting  from  the  simultaneous  detonation 
of  two  Comp  B  filled  steel  pipe  bombs  (15,2  cm  long,  4.83  cm  diameter,  3.2 
mm  thick  wall)  is  seen  in  the  end-on  flash  X-ray  photograph,  Fig.  1.  The 
details  of  the  jets  formed  by  this  mechanism  vary  considerably  from  shot  to 
shot  (largely  because  the  pipe  bombs  are  not  precision  items).  Typically, 
the  jets  travel  at  about  2  km/s  with  the  lowest  velocity  material  effective 
in  penetrating  metal  plates  traveling  at  about  0.5  km/s.  Material  in  the 
jet  is  confined  to  a  fan  of  material  about  0.1  radian  in  angular  width. 

The  design  of  an  optimum  suppressive  structure  requires  that  this  type 
of  jet  be  either  interrupted,  dissipated  or  inhibited.  We  describe  below 
the  model  used  to  calculate  the  effect  of  jet  collision  with  elements  of 
the  suppressive  structure. 

I'll.  ARTILLERY  SHELL  INFLUENCE  INITIATION  EFFECTS 

Impact  of  explosively  accelerated  metal  on  artillery  shells  can  produce 
influence  initiation  by  any  of  the  following  mechanisms: 

(1)  Shock  loading  of  the  explosive  arising  from  metal  on  metal  impact 
loading  to  the  incidence  of  a  strong  shock  propagating  into  the  explosive. 

(2)  Adiabatic  compression  or  plastic  flow  heating  of  the  explosive  (or 
entrapped  gas)  due  to  compressive  deformation  of  the  impacted  shell  or 

to  fragment  penetration. 

(3)  Plastic  flow  heating  of  explosive  injected  into  shell  fissures 
formed  by  metal  on  metal  impact. 

(4)  Direct  exposure  of  the  explosive  fill  to  hot  fragments  in  the  explo¬ 
sive  environment  due  to  shell  wall  failure  or  penetration. 

Those  mechanisms  entail  different  quantitative  criteria  to  determine  if 
initiation  will  be  produced  in  a  given  test  situation.  These  include: 

(1)  The  time  integral  of  the  shock  pressure  p  squared,  /p2 3dt .  This 
quantity  can  be  compared  with  critical  values  for  different  explosives . 

(2)  The  maximum  pressure  to  which  the  explosive  is  subjected  during 

shock  loading,  p  ,  and  the  rate  of  loading  dp/dt  of  the  explosive  or 
fnsx 

entrapped  gas. 

(3)  A  shell  deformation  characteristic  measure  A, 

A  =  J dt  dt 
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Figure  1.  Flash  X-ray  photograph  of  two  Comp  B  filled  steel  pipe  bombs 

from  end-on  view.  The  pipe  bombs  are  15.2  cm  long,  4.83  cm  in 
]  diameter  with  3.2  mm  thick  walls  and  separated  by  a  distance  of 

;  6.35  mm.  The  flash  X-ray  was  taken  50.1  us  after  detonation  of 

•  the  rounds.  The  vertical  stripe  in  the  X-ray  photograph  is 

caused  by  the  metal  jet,  which  of  course,  is  traveling  with  a 
r<  velocity  gradient. 


B 


i 


f; 
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where  a  is  the  effective  strength  of  the  shell  wall  material  (taking 
account  of  moment  arm  effects  in  the  loading),  and  p  is  the  density  of 
the  material . 

(4)  Depth  of  penetration  of  impacting  fragments  or  jets  into  the 
casing  of  the  artillery  shell. 

In  the  present  report  we  will  only  be  concerned  with  the  last  of  these,  the 
penetration  depth  into  the  wall,  as  from  the  experimental  data  it  appears  that 
this  is  the  most  critical  measure  in  determining  whether  or  not  metal  jetting 
will  lead  to  HE  initiation, 

IV.  DESCRIPTION  OF  THE  INFLUENCE  SUPPRESSIVE  STRUCTURE  PARAMETERS  FOR 
THE  ARTILLERY  SHELL  PACKAGING  AND  STORAGE  SHIELDS 

Experimental  tests3’1*  have  demonstrated  that  the  use  of  multiple  sheets 
of  material  (perpendicular  to  the  line  joining  the  center  of  the  detonated 
or  donor  shell  and  the  protected  or  acceptor  shell)  proves  to  be  effective 
in  suppressing  influence  initiation.  Plaster,  containing  water  of  hydration, 
has  appeared  to  be  a  particularly  good  material  for  use  in  these  experiments.*4 
In  addition  the  use  of  explosive  material  sandwiched  between  layers  of  inert 
material  may  prove  to  be  a  good  material  to  shield  artillery  shells  from 
metal  jets.  As  such,  a  typical  shield  design  may  have  as  many  as  six  layers, 
each  of  which  must  be  designed  to  be  of  optimum  thickness  and  placement. 

Figure  2  shows  the  general  configuration  to  be  treated  in  the  present  paper. 
Since  plaster  contains  water  of  hydration,  there  is  the  possibility  for  it 
to  behave  as  an  active  material  in  these  structures.  The  high  velocity 
impact  of  metal  from  an  exploding  shell  givos  rise  to  shocks  in  the  plaster, 
producing  shock  heating  of  the  material.  This  heating  can  liberate  water 
from  the  gypsum  in  the  form  of  steam,  driving  the  material  as  though  it 
had  exploded. 

In  the  following  treatment  of  the  suppressive  structure,  Fig,  2  will 
be  used  as  the  basis  for  the  physical  configuration  to  be  treated  by  the 
model.  For  the  most  part,  materials,  dimensions,  etc.  will  be  represented 
parametrically.  Subsequently  the  model  will  be  exercised  to  determine  those 
characteristics  yielding  the  best  results  under  various  constraints  (such  as 
total  weight  or  volume). 
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DONOR  ROUNDS 


Figure 


ACCEPTOR  ROUND 


2,  General  configuration  of  shields  treated  in  the  present  paper. 

The  donor  shells  are  assumed  placed  in  the  worst  case  position 
relative  to  the  acceptor  shell.  Two  shields  each  of  up  to 
three  layers  of  any  thickness  are  assumed  placed  at  any  position 
between  the  donor  and  acceptor  shells.  The  middle  layer  of  each 
shield  is  assumed  to  consist  of  an  inert  material,  a  semiactive 
material  such  as  plaster  with  water  of  hydration,  or  an  explosive 
material . 


V.  SHOCK  ACTIVATION  AND  JET  PENETRATION  OF 
HYDRATED  OR  EXPLOSIVE  SOLIDS 

Use  of  gypsum  in  the  form  of  plaster  panels  has  yielded  some  experimental 
data4 ,  suggesting  that  under  shock  loading  the  material  responds  in  an 
explosive  manner.  Sufficiently  strong  shock  can  give  rise  to  the  heating  of 
the  granular  solid  material  comprising  these  panels.  The  heat  acts  to 
calcine  the  plaster  giving  rise  to  a  high  temperature  steam  under  elevated 
pressures.  As  a  result  some  of  the  original  shock  energy  can  be  dissipated 
in  the  plaster.  The  subsequent  explosive  expansion  of  the  plaster,  however, 
may  give  rise  to  shocks  in  neighboring  materials.  Shocks  propagating  through 
water,  however,  have  not  exhibited  any  increased  ability  to  attenuate  shocks 
from  donor  rounds.  Thus,  there  may  exist  an  effect  due  to  shock  induced 
calcination  that  involves  the  "solid  residuum"  serving  as  the  source  of 
the  heat  available  for  the  generation  of  the  steam. 

Let  us  assume  that  we  have  a  completely  inelastic  collision  between  a 
plate  of  areal  density  (volumetric  density  p^)  traveling  at  an  initial 

velocity  v1  and  a  stationary  plaster  plate  of  areal  density  a2  (volumetric 

density  p2)  containing  a  fraction  g  of  moisture  as  shown  in  Fig.  3.  After 

the  collision  the  velocity  (of  the  center  of  mass)  is  v2*  The  energy 

absorbed  per  unit  area  E  is  given  by 

E  -  Jj  [P21/cl  -  ?22/(o1  -o2)]  (1) 

where 

Px  •o1v1  (2) 

and 

P2  "  (*i  +  a2)v2  (3) 

The  specific  energy  e2,  per  unit  mass  in  the  second  (plaster’)  plate  before 

expansion  is 

e2  =  E/o2  (4) 

This  energy  gives  rise  to  the  heating  of  the  "solid  residuum"  of  specific 
heat  C  (which  for  plaster  has  a  value  of  0.210  cal/g°C),  and  the  water 

5 

of  specific  heat  C  (*0.48  cal/g°K)  that  is  liberated  as  the  gypsum  is 

vh2o 

heated.  The  heat  of  liberation  is  represented  by  h^.  For  gypsum,  hf  is 
58.1  cal/g;  for  water,  79,7  cal/g.  We  obtain 

e.  =  (1  -  g)C  AT  +  C  AT  +  ghf  (5) 

*  'h20 
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Schematic  of  an  inelastic  collision  between  an  inert  plate,  A, 
traveling  with  velocity  v^' impacting  on  a  plate  containing  water 

(for  example,  a  plaster  panel  containing  water  of  hydration)  con¬ 
tent,  B.  The  resultant  combination  of  plates,  C,  travels  away 
with  a  velocity  v2>  the  kinetic  energy  absorbed  by  the  inelastic 

collision  being  converted  to  heat  with  the  evolution  of  steam. 


(This  expression  represents  the  unexpanded  condition  of  the  heated  plaster 
so  that  the  heat  of  vaporization  does  not  occur  here.)  where  AT  *  T  -  T^ ,  T^ 

being  the  initial  temperature  of  the  plaster  and  T  the  temperature  to  which 
an  amount  e7  of  energy  heats  the  plaster.  Of  this,  an  energy  (specific) 

0  will  be  available  as  the  explosively  released  energy  of  the  steam.  Where 
gC  AT  is  the  energy  that  goes  into  heating  the  liberated  water,  the 

vh,o 

available  energy  is 


sc. 


'H,0 


ATn 


(6) 


where  n  is  the  work  efficiency.  Ke  will  use  for  n  the  ideal  thermodynamic 
efficiency  (for  a  Carnot  cycle)  although  the  actual  efficiency  is  lower.  The 
ideal  efficiency  is 


n 


(7) 


Where  T  is  high,  Eq .  (S)  serves  as  a  good  approximation.  Thus,  we  can  write 


-eH  J n  *  [a1  vj  -  (o1  +  o,)v^)/o,  -  2 [ (1  -  g)C  AT  +  ghJ  (8) 

2  “  s' 

or  substituting  for  AT  from  Eq.  (8)  and  solving  for  q  we  have 

*  cv 

2eHnO  *  n{J2talVl  "  (ol  +  a2)v2^  "  28hf}/U  +  (9) 

2  Vh2° 

The  quantity  2e^  Q  can  be  used  as  the  energy  constant  in  the  Gurney  equation 

in  order  to  calculate  the  subsequent  expansion  of  the  plaster  plate  as  it 
drives  adjoining  plates.  It  should  be  noted  that  the  assumption  of  complete 
inelastic  collision  leads  to  Eq.  (5)  giving  the  energy  available  for 
calcination  of  the  gypsum.  In  the  case  where  a  given  shock  propagates  into  a 
panel,  it  is  necessary  to  detei'raine  the  rate  at  which  shock  energy  is 
converted  to  heat  which  will  be  a  function  of  (1  -  g) .  The  value  for  dry 
gypsum  (water  of  hydration)  is  0.209.  Measurements  by  Gould  Gibbons,  Jr. 

(TBD  -  EEB)  yielded  a  value  of  0.290  for  wet  plaster  as  tested. 

Let  us  now  consider  what  occurs  when  a  high  velocity  jet  of  material 
impacts  a  material  such  as  plaster  which  explosively  releases  energy.  Figure 
-l  shows  the  flow  of  jet  material  and  target  material  in  the  coordinate  system 
moving  with  the  stagnation  point  in  a  target  under  steady  state  flow 
conditions.  Letting  v.  be  the  velocity  of  the  jet  and  V  the  velocity  of 

the  stagnation  point,  in  the  frame  of  reference  moving  with  the  stagnation 

point,  the  jet  approaches  the  stagnation  point  S  with  a  velocity  v,  -  V  while 

j  s 

the  target  material  approaches  with  a  velocity  V  .  To  maintain  steady  state 
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flow,  the  forces  acting  on  the  stagnation  point  must  be  balanced.  Therefore,  we 
can  write,  noting  that  the  shocked  plaster  gives  rise  to  a  pressure  p^  Q, 

-  vs)2  •  ’kW  *  PHjO  <>0) 


We  note  that  p^  Q  is  simply  equal  to  the  energy  density  of  the  steam,  q. 


as  given  by  F,q.  (9).  By  reference  to  Fig.  3  and  recognizing  that  in  jet 
penetration  to  an  incremental  depth  dP  an  incremental  length  of  penetrator 
dL  impacts  with  the  target  and  is  consumed,  we  see  that  the  quantities  o. 

and  o2  of  Eq. 

(9)  take  the  values 

Oj  -+  dLpj 

(11) 

and 

°2  dPpp 

(12) 

Where 

V1  “  Vj 

(13) 

we  have 

°1V1  "  C°l  +  °2)V2 

(14) 

or 

o,  dPp 

v2  ■  y»  ■  vja  ♦  sesj1 

(15) 

Therefore,  we 

have  the  integral  expression  for  P 

r  L 

P  ■ J  dL^p/Cl  +  a) 

(16) 

where 

0 

P  -  Pj/Pp 

(17) 

“  "  2Ph,0/pPVs 


(18) 


Where  p  and  a  are  constant,  we  can  write  (16)  as 


P  ■  L^p/(1  +  a) . 


Substituting  from  (13)  and  (14)  into  (9)  gives 

V  ■  >!,’k  (  ' 


2  °1 

°2  "  °l  +  a2  * 


where 


(1  -  R)C, 


k  -  [1  + 


8C, 


i]"1. 


HjO 


Defining 


we  can  write  (21)  as 


dL 


f  -  Ql/02  -  pjjp  , 


eH20  •  ^[Vj(~fvj  h  -  2ghf] 

Substituting  (24)  into  (19)  and  then  into  (18)  gives 

2,f2  +  f  -  1 


(20) 


(21) 


(22) 


(23) 


(24) 


(nMv  (--f‘ . y  -0  -6  ] 

(o  if  nk[v2(^-^-;-i-)  . 


(25) 
6  ]  <  0 


where  v 


v ,/V  and 
J  3 


«  -  2ghf/V2. 


If  we  take  the  effects  of  the  steam  evolution  to  bo  small  then 


(26) 


f  ■  3P7JT  Vp  *'/p  '  t27) 

With  this  value  as  an  approximation,  Eq.  (25)  can  be  solved  iteratively 
allowing  the  integration  of  Eq.  (16).  It  may  bo  seen  that  Eqs.  (20)  and 
(28)  abridge  the  usual  penetration  density  law  for  hypervelocity  target 
penetration,  giving  rise  to  a  term  that  increases  the  effective  density  of 
the  target  material.  Thus,  in  the  case  of  impact  activated  materials  the 
usual  hypervelocity  penetration  density  law  differs  from  that  appropriate  for 
inert  materials. 
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VI.  DESIGN  OPTIMIZATION  FORMULATION  FOR  INFLUENCE  ATTENUATION 
PANELS  TO  SUPPRESS  SHELL  -  SHELL  IMPACT  JET  INDUCED  INITIATION 

Assume  a  jet  formed  by  the  impact  of  two  simultaneously  initiated  shells, 
as  shown  in  Fig.  1,  attacks  an  influence  suppression  structure  of  the  kind 
shown  in  Fig.  2.  Figure  5  illustrates  the  specific  jet  and  suppressive 
structure  concerned  together  with  parameters  representing  the  important 
dimensions  and  quantities.  Assume  the  jet  formed  by  the  impact  of  the  shells 
has  an  initial  length  £  of  density  pQ,  the  density  of  the  shell  walls,  and 

that  subsequently  the  jet  increases  in  length  linearly  with  the  difference 
between  the  jet  tip  velocity  ,rt  and  jet  back  velocity  v^.  We  assume  also 

that  the  jet  width  increases  from  an  initial  value  e  with  a  velocity  av(x) 
where  a  is  a  constant  equal  to  the  angle  of  spread  in  jet  material  measured 
in  radians  and  v(x)  is  the  forward  velocity  of  the  jot  at  any  distance  x 
along  the  initial  jet  length.  The  density  of  the  jet  is  assumed  to  decrease 
inversely  with  the  product  of  the  fractional  increase  in  width  and  local 
stretching  in  length. 

We  assume  a  linear  velocity  gradient  given  at  any  point  x  by 

v  “  vb  +  (vt  "  Vb)x/A  (28) 

where  l  is  the  length  of  the  jet  at  which  tho  density  i3  that  of  the  steel. 
The  time  for  material  to  reach  any  distance  L  will  be  referred  to  the  time  of 
initial  jetting  (virtual  origin  of  the  jet  at  x»0) .  Thus, 

t (x)  -  L/v(x).  (29) 

The  jet  element  of  length  dx  at  x  increases  in  length  to  dx'  as  of  its 
arrival  at  L  as  given  by 

dx'  ■  dx  +  tdv  (30) 

Differentiating  Eq.  (28)  to  substitute  for  dv  and  using  (29)  to  substitute 
for  t  in  (30)  gives  on  integrating 

x'  «  [1  +  CXL/  (1  +  C^x)  ]dx  (31) 

-  I  *  L  ln(l  +  CjJt) 

where 

C1  "  (vt  '  vb)/aV  (32) 

This  gives  the  length  of  the  jet  at  any  distance  L  from  its  origin. 
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Figure  5 


DONOR  ROUNDS 


ACCEPTOR  ROUND 


Schematic  of  tho  jet  formed  by  simultaneously  initiated  artillery 
shells  and  antifratricide  panels  of  the  type  shown  in  Fig,  2 
defining  the  important  parameters  governing  the  interaction.  Jet 
tip  and  base  velocities  vt,  vb,  width  e,  angular  spread  a,' 

initial  length  l  and  jet  density  p.,  panel  densities  p  ,  p.,  p„, 

, ]  R  D  C 

pfll,  PjjM  and  o  ,  as  well  as  the  spacial  configurational  parameters, 
c1,  h,  d,  dg  and  L  are  shown.  Density  of  the  target  is  pt . 
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The  penetration  of  the  unimpeded  jet  at  L  in  a  target  (such  ag  steel)  of 
density  must  be  calculated  in  terms  of  the  density  of  the  jet  as  it 

arrives  at  the  target.  Since  the  jet  initially  forms  a  sheet  of  metal,  its 
density  drops  as  it  expands  laterally,  i.e.,  normally  to  the  surface  of  that 
sheet,  and  as  it  lengthens.  Thus,  we  have  for  the  density: 


°j(K't)  'so  .  *  .v(») 


m  (air) 


L  C.L  V1 

"  p°  I1  +  l  +  CjT  j 


e  +  avt 


The  penetration  factor  is  y,pj/pt .  Thus,  without  influence  suppression 
panels,  the  penetration  P  in  the  target  at  L  will  be 


a  _ r  Cjd  +  pcxj) 

70  dxVpj/pt  1  *  — rrepr — 


where  we  have  included  the  penetration  into  the  target  L  +  P(x)  as  a 
correction  to  more  accurately  give  the  distance  to  the  target.  This 
correction  is  to  be  included  in  the  expression  for  p^  as  well,  but  only  for 

the  longitudinal  extension  of  the  jet: 


Pj(x.t)  ■  p( 


/  cta  y wA-1. 

I  T+UjX  I  e  +  al 


This  gives  for  the  penetration  into  the  target  in  the  absence  of  suppressive 
layers : 


PWpt/dx  1  + 


Ca  CL  +  Px) 


+  CjX 


where  the  notation  Px  means 


px  i  P(x)  .  *  ~T~7 *  t  l)1/2(37> 

We  now  calculate  the  effects  due  to  the  presence  of  two  suppressive 
structure  panels  through  which  the  jet  must  pass  before  encountering  the 
target  plate.  Assume  panel  1  consists  of  the  three  layers  of  thickness 
a,  c,  and  b  as  shown  in  Fig.  6  We  have  for  the  losses  of  jet  length  M  , 


Mc,  AAfa 


and 


1/2  (38) 


c  * 


, —  r  ci(h  +  pttn1/2  /  /  „  \ 

*  **■  ['  *  -TT-c— J  /  ('  *  I  h) 

£-aaa  ' 

_ _  ,-A-AA,,  r  c.  (h+a+PJ-]1^2 

Vp»/pi /  d*  [“  TTcjx  '  ?  (h  *  “’]  <39> 


A-AA  -AA 
a  c 


A-AA  -AA 

a  i 


b  “  *lp0/pb 


dx 


1  + 


A-AAa-AAc-AJ^ 


Cjfhtoto-tP^I1^ 

>'  *  *,*  J/L 


a  1 1/2 

1  +  2.  (h+a+c)  I  (40) 


where  Pft,  P  ,  and  P^  are  defined  as  the  instantaneous  penetration  as  in  Bq. 
(37).  The  density  is  from  (20),  the  "effective"  value 


pj  ■  pc(l  +  a)  (41) 

-  '*>  . 

For  high  velocity  jet  penetration,  stagnation  point  temperatures  are  high,  so 
that  the  efficiency  n  can  be  nearly  unity.  In  this  case  (41)  becomes 

which  is  to  be  solved  recursively.  The  total  loss  in  length  passing  through 
panel  1  is  therefore  AAj  where 


A  A,  ■  AA  +  AA  +  AA, 

1  a  c  b 


(42) 


During  the  jet's  penetration  of  panel  1,  a  shock  builds  up  in  the  Jet 
as  shown  in  Fig.  4.  When  the  jet  exits  panel  1,  the  rarefaction  wave  from 
the  newly  formed  free  surface  of  the  jet  propagates  into  this  region  of 
high  pressure,  disrupting  the  jet  material.  Let  us  obtain  an  approximate 
expression  for  the  loss  of  jet  material  due  to  this  effect.  There  are  two 
cases  to  be  considered.  First,  if  the  penetration  is  supersonic  an  amount 
proportional  (and  approximately  equal)  to  the  jet  thickness  is  lost,  i.e.,  and 
amount  ke^,  where  k  is  a  proportionality  constant  approximately  equal  to  unity. 
Since  this  should  be  expressed  so  that  a  zero  thickness  panel  removes  no  jot 
material  we  will  write  for  the  loss  length 
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(43) 


yy  «  kej/  [1  +  6jpj  / (apft  +  bpb  +  cp^] 

where  e1  Is  the  jet  width  given  in  terms  of  the  initial  width  by 

e1  *  e/ [1  +  “  (h  +  a  +  b  +  c) ]  (44) 

and  p.  is  the  density  of  the  jet  on  exit  from  the  first  panel.  In  this 
31 

length  y^,  shocks  of  pressure  p,  given  by 


plpj, 


Pi  - 


2M3!  +  %/Pj^2  1 


(45) 


where  should  be  taken  to  be  the  average  density  within  a  mass  column  equal 
to  e^,  act  on  the  mass  column,  In  essence  this  means  unless 

b<<o  or  pb<<pj.  Also,  Pj  should  be  taken  to  be  the  jet  density  at  exit  from 
panel  1;  * 

C^h+a+b+cjV1  e 

Pj.  "  ‘  "  ' 


po  I1  +  i'+C1'(i-Ai1)'l  eWfa-AijV 


(46) 


Here  Vj  is  the  corresponding  velocity  at  x  “  A-  AAj  as  given  by  Bq.  (38) . 
Now  the  pressure  p^  causes  the  jet  to  spread  laterally  and  thus  drop  in 
effectiveness  against  the  next  plate.  This  loss  in  effective  length  Ay1 
is  calculated  from  density  considerations  to  be 


•*i  [1  -  tJTT  +  ad/e)/  (1  +  afjd/ef  ] 


(47) 


whore 


ov, 


1  ovj  +  (p1/p^)1^5' 


(48) 


In  tho  case  of  subsonid  penetration  of  the  punel,  the  shock  zone  in  the  jet 
will  increase  with  penetration  distance.  However,  the  strength  of  the 
shock  will  decrease.  Thus,  Bq.  (47)  will  serve  as  an  approximation  of  the 
jot  length  loss  in  that  case  as  well.  This  length  loss  is  equivalent  to  an 
initial  length  decrement  AAj'  given  by 
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f  JC-flJt  J 

1  "I 

*'»_(■  AliAO  » 


dx  |1  + 


ni- (aji+aHj  » D 


C.  (h  +  a  +  b  +  c  +  d) 

1+  -i - - - - - 

1  +  CjX 


At  panel  2  the  length  loss  is  to  be  computed  in  the  same  fashion  as  for 
the  first  panel.  We  have  in  the  appropriate  order: 


b>  ■  \VPI 


r  c  (L  p  ,-.1/2  / 

j  *[ltJkrhr-\  /  "i1*] 

A£1+Ai1’)-Ailb,L  J  /  J 


where 


L^^h  +  a+  b  +  c  +  d 


_  f  r  C.(L+bUP  .)■  '  /r 

yfc/p2j  dx  jl  + - JTTCjX -  /  1  +  e 


(L2+b')  (52) 


V^a 


«,-(AJ,1+AS.1')-A«.b,-Alc,  -1/2/ 

-/*  r  C,(L,+b'+c'+Po,)‘  V 

•/  dx  [1+  1  f v j/[i  * 5- (vb,+c,l 


where 


Li2  “  Aila'  +  A4b'  +  &2,c* 


As  in  the  case  for  Eq.  (41),  we  write  for  p2 


P2  "  Pci  (1  +  «')  -  pc,  <  1  +k  [v 


rrr^  - 4  lj 


Again  we  calculate  the  effective  loss  in  jet  length  due  to  shocks  in  the  jet 
on  exit  from  the  second  panel.  Corresponding  to  Eq.  (43)  to  (49)  we  have 


y2  -  ke2/[l  +  e5p^/(a'pn,  +  b'pM  +  c'p,)] 


2Mj2  Ma 
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with 


&2  =  e/[l  +  e  (L>2  +  a'  +  b'  +  c')]. 


(57) 


The  shock  pressure  p2  is 


P2 


2(  y/  C>2  +  /Pi  ) 
"  J  2 


2  V2 


(58) 


where 


-1 


pj2  "  pofl  + 


CjCLj+a'+b'+c*) 

1 +C  x  ( £  -  A  &  x  >  A  Xp*  -A^2V  I  e+a(£-AJlY-AJ!,1,-A2,2) 


(59) 


and  v2  is  the  velocity  at  x  *  1  -A)l2  as  given  by  Eq.  (38).  The  spread 

in  the  jet  produces  the  effective  loss 


Ay2  *  -M1  -  ^  +  ah' /a) /Cl  +  af'h'/e)  '  ] 


(60) 


where 


av. 


f' 


av2  +  (P2/pj2} 


172 


(61) 


Therefore,  A£2',  the  effective  loss  in  initial  rod  length  due  to  the  gap  h' 


is 


.1- (AH, +Ad' ) -A£-  T  C. (L^+a' +b' +c' +h' )" 

Ay2  =  ^  dx  ^1  +  - j  (62) 

J  A-fA^+AJip-Atj+A^) 

Finally,  the  remaining  rod  material  penetrates  into  the  target  a  distance 

il-  (AS^+AR^  1 )  -  (AS,2+A1.2 ') 


P'  ■  w  — 


dx 


C1<L+PV 

1  +  Cxx 


1  + 


1/2 

(63) 
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VII.  COMPUTATIONS 


A  computor  code  designated  "Main"  to  compute  P'  has  been  developed  and  is 
given  in  Appendix  A.  This  code  has  been  used  to  compute  shield  effectiveness 
curves.  The  effectiveness  of  a  single  homogeneous  panel  to  suppress  jet 
induced  influence  is  shown  in  Figs.  6-9.  The  nondimensional  penetration  (in 
units  of  donor  round  radius)  P'  as  a  function  of  the  nondimensional  panel 
thickness  (also  in  units  of  donor  round  radius)  is  plotted  for  the  case  of 
separation  distance  dg  equal  to  1.0  R  (Fig.  6),  2.0  R  (Fig.  7),  3.0  R  (Fig.  8), 

and  4.0  R  (Fig.  9).  The  shield  panel  is  taken  to  be  equidistant  between  the 
donor  and  acceptor  rounds.  Panel  density  pc  is  taken  as  a  parameter  for  Figs. 

6-9.  Values  of  specific  parameters  have  been  taken  to  provide  an  optimum 
fit  to  experimental  data  obtained  from  pipe  bomb  tests  as  discussed  below.4 
X-rays  made  during  these  pipe  bomb  tests  were  employed  to  obtain  the  jet 
characteristic  parameters'  values.  Specific  values  used  are: 

Jet  Characteristics:  pQ  ■  8.0  g/cm3 

l  =  0.2 
e  =  0  .4 
a  =  0.1 
vt  =  2.0  km/s 

vfe  *  0.5  km/s 

Target:  Pt  *  8.0  g/cm3 

Other:  k  «  2.0 

vs  «  3.0  km/s 
g  =0.0  and  0.209 

The  results  obtained  for  the  two  values  of  g  were  neatly  the  same.  The 

values  of  jet  characteristic  parameters  are  dependent  on  the  donor  rounds 

used.  For  rounds  that  do  not  scale  to  the  approximate  values  of  the  pipe 
bombs  used  in  Ref.  4,  scaling  test  should  be  made  before  the  results  of 
Figs.  6-9  are  used  for  applications.  The  curves  offer  no  surprises,  but 
simply  show  that  increased  thickness  and  increased  density  both  decrease 
monotonically  the  penetration  in  an  acceptor  round.  The  curves  can  be  used 
to  obtain  the  requirements  to  just  stop  a  jet  from  penetrating  into  the 
acceptor. 


PVR 


DONOR  ROUNDS 


Figure  6.  Nondimensional  penetration  P'/R  versus  nondimensional  shield 

thickness  C/R,  where  shield  is  placed  at  midpoint  between  donor 
and  acceptor  rounds  (or  witness  plate,  as  shown  in  inset)  separated 
by  1.0  R.  Shield  density  p  is  parameter  for  the  family  of  curves. 


PVR 


DONOR  ROUNDS 


C/R 


Figure  7.  Nondimensional  penetration  P'/R  versus  nondimensional  shield 

thickness  C/R,  where  shield  is  placed  at  midpoint  between  donor 
and  acceptor  rounds  separated  by  2.0  R.  Shield  density  pc  is  a 
parameter  for  the  family  of  curves. 


PVR 


DONOR  ROUNDS 


C/R 


Figure  8.  Nondimens ional  penetration  P'/R  versus  nondimensional  shield  thick¬ 
ness  C./R,  where  shield  is  placed  at  midpoint  between  donor  and 
acceptor  rounds  separated  by  3.0  R.  Shield  density  pc  is  parameter 
for  the  family  of  curves. 


PVR 


DONOR  ROUNDS 


Figure  9.  Nondimens ional  penetration  P'/R  versus  nondimensional  shield 

thickness  C/R,  where  shield  is  placed  at  midpoint  between  donor 
and  acceptor  rounds  separated  by  4,0  R.  Shield  density  p  is 
parameter  for  the  family  of  curves,  c 


Figure  10  shows  the  results  obtained  for  the  nondimens ional  penetration 
depth  as  a  function  of  the  separation  distance  D/R  for  two  panels  placed 
symmetrically  about  the  midpoint  between  the  donor  and  acceptor  (or  target) 
rounds.  The  result  is  obtained  for  several  densities  and  for  panels  0.53  R 
thick  (this  corresponds  to  panels  1/2"  or  1.27  cm  thick  in  the  pipe  bomb 
experiments) .  Separation  between  donor  and  acceptor  rounds  Is  taken  to  be 
4.0  R  in  these  calculations.  These  curves  show  that  the  penetration  versus 
separation  distance  drops  as  the  separation  increases.  This  is  easily  under¬ 
stood  as  an  effect  of  jet  elongation  with  travel  distance.  As  a  result,  it 
is  concluded  that  for  a  fixed  amount  of  buffer  material  it  is  better  to  separ¬ 
ate  the  material  into  two  panels  placed  close  to  or  against  the  shells.  Fig¬ 
ures  11-14  incorporate  the  optimum  shield  material  distribution,  i.e.,  separ¬ 
ation  of  the  shield  material,  into  two  panels  as  indicated  by  the  results  shown 
in  Fig.  10.  Comparison  with  Figs.  6-9  shows  the  improvement  afforded  by  such 
an  arrangement  of  shield  material.  Note  that  the  curves  in  Figs.  6-9  give  the 
same  value  at  zero  total  shield  thickness  and  for  shields  completely  filling 
the  donor-acceptor  round  gap  as  the  corresponding  curves  in  Figs.  11-14.  At 
points  in  between,  however,  the  curves  in  Figs.  11-14  dip  below  their  counter¬ 
parts  In  Figs.  6-9.  These  results  show  that  the  improvement  noted  in  Fig.  10 
as  to  the  optimum  placement  of  shield  material  under  the  symmetry  constraint 
is  general. 

Figure  15  gives  data  showing  that  a  further  improvement  in  buffer  material 
efficiency  is  achieved  using  a  "foam"  (a  low  density  material)  filling  the  gap 
between  the  shells.  The  curve  gives  the  penetration  as  a  function  of  the  den¬ 
sity  of  the  "foam"  buffer  filling  the  entire  gap.  It  is  seen  that  the  pene¬ 
tration  in  the  target  plate  is  zero  at  a  density  of  about  0.21  g/cm3  for  a  gap 
of  4. OR.  This  is  equivalent  to  a  1.05  R  thick  panel  of  density  O.BO  g/cm3 
which  placed  at  the  midpoint  will  not  protect  the  target  plate  (penetration 
P'/R»0.15);  divided  into  two  0.525  R  thick  panels  the  penetration  is  P'/R-0.13. 
Therefore,  the  best  panel  design  appears  to  be  a  low  density  "foam"  filling 
the  available  gap, 


VIII.  COMPARISON  WITH  EXPERIMENTAL  RESULTS 

An  extensive  series  of  shield  tests  using  4.82  cm  (1.9")  diameter  pipe 
bombs  has  been  run  employing  a  wide  variety  of  shields.  The  data  obtained 
in  those  tests  has  been  used  to  validate  the  "Main"  code.  Figure  16  shows 
a  comparison  between  the  experimental  and  theoretical  code  results.  In  Fig. 

16  the  penetration  depth  (maximum  depth  in  the  target  plate)  is  used  for  the 
ordinate  while  the  Individual  test  code  numbers  (these  are  the  samo  code 
numbers  used  to  designate  the  tests  in  Ref.  4)  are  listed  along  the  abscissa. 
The  order  of  the  tests  listed  along  the  abscissa  is  not  important.  As  listed, 
however,  they  imply  a  decreasing  function  of  penetration  with  increasing 
shield  material.  Boneuth  each  entry  is  a  small  picture  indicating  the  test 
arrangement  of  the  shield  along  with  the  shield  thickness  in  mm  together  with 
symbols  indicating  shield  materials  used.  Additional  test  details  may  be 
found  in  Ref.  4.  All  tests  employed  Comp  B  filled  4.82  cm  diameter  pipe  bombs 
15  cm  in  length  placed  4  radii  above  a  steel  witness  plate.  The  reference 
case  shows  the  results  for  no  shield  between  the  two  pipe  bombs  and  the 
witness  plute,  The  error  bars  (obtained  from  the  tost  data  for  three  shots) 
indicate  the  considerable  variability  of  the  results  obtained  in  these  tests. 
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PVR 


DONOR  ROUNDS 


Figure  10.  Nondimen sional  penetration  P'/R  versus  nondimensional  panel 
separation  distance  D/R  for  two  symmetrically  placed  0.53  R 
thick  panels  in  a  4.0  R  space  between  donors  and  target  plate  for 
various  panel  densities.  The  inset  shows  the  geometry.  Those 
results  indicate  that  if  symmetry  must  be  maintained  in  packaging, 
the  best  placement  of  shield  material  is  midway  between  donor  and 
acceptor  rounds. 
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Figure  11.  Nondimensional  penetration  P'/R  versus  nondimensional  total 

shield  thickness  (C  +  C')/R,  where  the  shield  consists  of  two 
panels  symmetrically  placed  about  the  midpoint,  one  against  the 
'donor  rounds,  one  against  the  acceptor  rounds  (or  witness  plate, 
as  shown  in  the  inset).  Donor-acceptor  separation  is  1.0  R. 
Shield  density  p„  is  parameter  for  curves, 
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Figure  12.  Nondimens ional  penetration  P'/R  versus  nondimensional  total  shield 
thickness  (C  +  C')/R,  where  the  shield  consists  of  two  panels 
symmetrically  placed  about  the  midpoint,  one  against  the  donor 
rounds,  one  against  the  acceptor  rounds  (or  witness  plate  as 
shown  in  inset).  Dono- -acceptor  separation  is  2.0  R.  Shield 
density  p  is  parameter  for  curves. 
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Figure  13.  Nondimen sional  penetration  P 1 /R  versus  nondimens  Iona  1  total  shield 
thickness  (C  +  C')/R,  where  the  shield  consists  of  two  panels 
symmetrically  placed  about  the  midpoint,  one  against  the  donor 
rounds,  one  against  the  acceptor  rounds  (or  witness  plate,  as 
shown  in  inset).  Donor-acceptor  separation  is  3,0  R.  Shield 
density  p„  is  parameter  for  curves, 
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Figure  14.  Nondimensional  penetration  P'/R  versus  nondimen sional  total  shield 
thickness  (G  +  C')/R,  where  the  shield  consists  of  two  panels 
symmetrically  placed  about  the  midpoint,  one  against  the  donor 
rounds,  one  against  the  acceptor  rounds  (or  witness  plate  as 
shown  in  inset).  Donor-acceptor  separation  is  4.0  R.  Shield 
density  p„  is  parameter  for  curves. 
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Figure  15.  Nondimonsional  penetration  versus  density  of  shield  material  for 
a  "foam"  filling  the  space  between  the  donor  rounds  and  target 
plate  (acceptor  rounds)  for  a  gap  of  4,0  R, 
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JET  CHARACTERISTICS  :  Zo=0.5cm 
VELOCITY  BACK/TIP  :  0.5-2J0km/s 
WIDTH/SPREAD  :  1.0cm.  0.1  rod 


In  tests  26T  and  30B  Detaaheet  was  employed  between  thin  sheets  of  steel 
to  form  the  shield.  The  code  computed  shielding  assuming  both  detonation  of 
the  Detasheet  and  no  detonation.  The  results  assuming  no  detonation  are  seen 
to  agree  satisfactorily  with  the  observed  results.  At  the  right  side  of  Fig. 

16  are  located  all  the  tests  in  which  a  shield  panel  is  placed  against  the 
donor  rounds.  It  will  be  noticed  that  these  six  tests  yielded  the  best 
shielding  (least  penetration  into  the  witness  plate)  results.  This  is  in 
general  agreement  with  the  predictions  of  the  code  as  indicated  by  the.  results 
in  Fig.  10.  Overall  the  theoretical  results  given  by  the  code  correspond  well 
with  the  experimental  results  in  Fig.  16.  (The  reader  should  keep  in  mind 
that  in  some  cases  the  penetration  data  is  a  measurement  taken  for  two  or 
three  fragment  craters  in  the  witness  plate  that  may  be  no  more  than  2-3  mm  in 
diameter.)  The  experimental  data  must  be  regarded  as  having  significant  error 
bars  reflecting  the  variability  of  these  test  results.  Within  these  limita¬ 
tions,  therefore,  the  model  works  well. 

The  results  obtained  in  matching  the  experimental  and  theoretical  data 
indicated  that  the  effects  observed  were  due  to  the  mass  column  of  shield 
material,  and  placement  of  that  material  alone.  The  code  indicated  no  signi¬ 
ficant  effect  due  to  thu  use  of  hydrates  in  the  Bhield.  Several  shots  to  test 
this  hypothesis  have  been  carried  out  by  G.  Gibbons4.  In  these  shots  equal 
mass  columns  of  steel,  aluminum,  plaster,  and  water  were  used  as  shields. 

Damage  to  the  witness  plates  was  almost  identical  in  these  four  shots,  although 
in  the  shot  using  steel  for  the  shield,  the  shield  material  itself  produced 
some  pitting  in  the  witness  plate  (these  pitB  were  located  to  either  side  of 
the  area  usually  damaged  by  the  jet  from  the  donor  rounds) . 

IX.  CONCLUSIONS 

If  two  neighboring  artillery  shells  in  an  ammunition  pallet  are  simul¬ 
taneously  detonated i a  high  velocity  metal  jet  is  formed  that  can  detonate 
other  neighboring  shells,  even  if  these  shells  are  protected  by  a  significant 
thickness  of  armor  (for  example  15  cm  of  aimor  can  be  penetrated  by  the  Jet 
from  two  simultaneously  detonated  155  mm  HE  artillery  rounds).  Equations  to 
model  the  interaction  of  Bhield  panels  with  such  jets  have  been  derived  and 
employed  to  write  a  computer  code  to  calculate  penetration  depth  of  the  jet 
after  passing  through  up  to  two  panels  consisting  of  sandwiches  of  up  to  three 
layers  of  material.  In  the  code  one  layer  in  each  panel  (the  middle  layer  in 
three  layer  panels)  can  consist  of  an  explosive  or  a  hydrated  material  that 
may  be  activated  to  yield  enhanced  pressures  under  the  shock  loading  of  the 
jet. 


Computations  using  the  code  have  indicated: 

(1)  The  primary  effect  produced  by  inert  or  hydrated  materials  depends  on 
the  mass  column  of  material  in  the  panel  and  panel  placement.  Water  of  hy¬ 
dration  produces  only  minor  effects. 

(2)  The  optimum  placement  for  the  shield  panel  to  suppress  jets  is 
against  the  donor  Bhells. 
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(3)  Under  symmetry  requirements  (i.e.,  where  all  shells  oust  be  equally 
shielded)  the  best  placement  employs  two  shields  of  equal  thickness,  one 
against  the  donor,  one  against  the  acceptor  rounds. 

(4)  Use  of  explosive  material  sandwiched  in  the  shields  can  significantly 
enhance  the  effectiveness  of  the  shields.  (Such  material,  however,  may 
enhance  impact  shock  initiation,) 

(5)  The  best  symmetric  shield  design  (minimum  weight  basis)  for  inhibiting 
fratricide  jets  would  employ  a  single  panel  filling  the  entire  available 
space  between  donor  and  acceptor  rounds  with  the  lowest  density  material 
capable  of  inhibiting  the  jet  (as  given  by  the  curves  in  Figs.  6-9)  in  the 
available  space. 
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PROGRAM  MA IN ( INPUT* OUTPUT , TAPE  5- INPUT.TAPE6- OUTPUT*  TAPE9 ) 

DIMENSION  GG ( 20) ,  DD(20>,  RRH0(20) 

DIMENSION  L  A  B  ( A )-»  PL C <  50 >*  PLPP(50) 

C 

C  ***  FORMATS 
C 

1  FORMAT!*'  GAP  -  ■  *  Ell  .4# /> 

2  FORMATt 3X* 'ICC  ■  ', I 10»4X, * PLC  ■  '  • Ell .4* 3X, • PL PP  «'*E11.4,//> 

A  FORMAT ( A10  ) 

5  FORMATt A10,10X,* LL  »',I3) 

10  FORMAT (1H1 ) 

15  FORMATt 1H0) 

20  FORMAT (3E10.3) 

25  FORMAT (IS) 

30  FORMAT  (  AX,  'A  -* *  E 11 .  A,  6X,  •  B  ■ 1  ,£11.4*5X»  •  BP  ,E11.  A,  5X, '  AP  -«,E11 
1.4»6X*'C  «*,E11.4»1X»'CP  -  •,Ell.A*/#AX,'t)  ■",E11.A#»  RADIUS  •»»E11 
2.  A,6X» 1  H  -*,£11 <  4>5X ,  *  HP  ■ '»  Ell .4,5X, ' AL  , E11.4, 3X* • RHOZ  -',E11. 
3A» /»  AX, " E  ■ ' *  Ell . 4* 5X* • VT  E11.4,5X, » VB  ■  •  , E 11. A* 2X* • AL PHA  «i*El 
A'l  .A, 3X*  •  RHOT  Ell.  4<*3X,»  RHOA  Ell. A, /,1X,"RH0B  ■ **€11 . 4*  3X," RH 

50C  *'*E11.4»2X,'RH0AP  ■« , E 11 . 4* 2X, ' RHOB P  «» , Ell . 4,2X* • RHOC P  -«,E11 
6 .  A,  6X, '  G  ■  ' »  E 11 .  A, /,  3X,  '  VS  ■ '  ,E11  .A#  5X*  '  AK  ■  ».,  Ell. A,  5X**HF  -  • » E 1 1. 
74 ,/) 

AO  FORMAT ( 9X, 'P  ■• ,  Ell .4* 5X* • Cl  ■ ' , E 11. A, 3X, • BETA  ■»,E11.A) 

10  FORMAT ( 1H  ,*  CAPL  ■» ,ill.4,5X» *  PA  ■' »E11. 4, BX*  * PB  ■ • , Ell . 4*5X, • PC 
1- ' » Ell . 4*5  X, *  VI  , Ell .A, 3 X, '  RHOP  -SEll.A*/,*  RH01  -**£11 .A* 2X, * 

2RH0J1  ■*, E 11 • A, 3 X," PRES  -',E11.4) 

60  FORMAT  (  3X»  •  F  ■'•*  E 11 .  A,  5X, »  Y1  • ' * E 11 .4, 2X» » DEL Y1  ■ ' ,  E 11 .  A,  3 X,  *  P2 
1»  Ell  .A, *  DELL1P  ■  ' , Ell .A ) 

62  FORMATt*  RH02  »*»E11.A) 

70  FORMAT  (  3X,  *PBP  •  • , E 1 1. 4,  5X , 'B P  »*,E11.4,'  DELL  BP  -',E11.4) 

80  FORMATt AX* «Y2  • ' » Ell . A* 5X» • V2  ■*• Ell . A, 3X, • PRES  - '* Ell . A* 6X,«F  »«, 
lEll. A, 2X>*  DELY2  Ell .4,5X, ' P4  «'*E11.4»/*'  DELL2P  -**E11 .4*5X, 'F 
2P  ■' ,E11»A ) 

B2  FORMAT  <  1H  , "DELLA  -Sill. 4,'  DEL  LB  -*,E11.A**  DELLC  -**E11.4) 

90  FORMATt AX, 'PP  ■' *  Ell . A, 2X, ' DELL1  «*,E11.4*'  DELL1P  ■ *,€ 11 . 4* 2X, • DE 
1 LL2  ■*',  Ell  •  A*  *  D1LL2P  ■«  ,E  11.  A*  3X*  « PPIN  ■•»E11.A) 

92  FORMAT 1 1H  ,’DELLl  «»,E11.A) 

9A  FORMATt  3X,  «PCP  -•  ,E11.  A,  5X, ‘CP  «*',E11.4*«  DELLCP  -*,E11.4) 

96  F0RMAT13X*  'PAP  •" , Ell . A, 3X," AP  ■•,E11.4*'  DELLAP  ■* ,E11 .A, 2X, ' DELL 
12  »',E11.A) 

98  FORMATt*  RHOPP  «'* ,  El  1 .  A,  2X  ,  « RHO  J  2  ■  *  ,  Ell .  A,  5X,  *  V2  ■»,E11.4) 

C 

C  ***  CALCULATIONS 
C 

CVH20  ■  2.007E7  $  CV  «  0.78E6  $  VDET  «  7.0E5 
WRITE! 6* 10) 

BX  ■  1.0  S  BY  ■  1.0 

LAB(l)  -  10H  TONTHMA  $  LAB (2 )  ■  10H  X9742**309 
LAB(3)  -  10H  PLOT  S  LABIA)  ■  10H  SCOOP 
XPAGE  ■  8.0  t  YP AGE  •  8.0 
CALL  PLTBEG(XPAGE,YPAGE,1.0,9,LAB) 

CALL  PL OT ( BX* BY* -3 ) 

BX  -  0.0  S  BY  -  0.0 
CALL  FACT0R(0.3) 

C  DO  3000  IC  -  1,2 
DO  3000  IC  ■  2,2 

READ (5*  20)  A*  B,  BP*  AP,  CP,  RADIUS*  HP,  AL,  RHOZ,  E*  VT*  VB, 

1  ALPHA,  RHOT*  VS*  AK,  HF 

A  -  0.0  *  B  -  0.0  $  AP  -  0.0  $  BP  ■  0.0  *  CP  •  0.0  $  HP  -  0.0 
RADIUS  >  1.0  $  AL  ■  0.20  Si-  O.AO 
READ (5,25)  ICOUNT 
R EAD (  5>,  A )  TITLE 

DO  700  LL  -  1, ICOUNT  43 


MAIN 

2 

MAIN 

3 

MAIN 

A 

MAIN 

5 

MAIN 

6 

MAIN 

7 

MAIN 

8 

MAIN 

9 

MAIN 

10 

MAIN 

11 

MAIN 

12 

MAIN 

13 

MAIN 

1A 

MAIN 

15 

MAIN 

16 

MAIN 

17 

MAIN 

18 

MAIN 

19 

MAIN 

20 

MAIN 

21 

MAIN 

22 

MAIN 

23 

MAIN 

24 

MAIN 

29 

MAIN 

26 

MAIN 

27 

MAIN 

28 

MAIN 

29 

MAIN 

30 

MAIN 

31 

MAIN 

32 

MAIN 

33 

MAIN 

3A 

MAIN 

35 

MAIN 

36 

MAIN 

37 

MAIN 

36 

MAIN 

39 

MAIN 

AO 

MAIN 

A1 

MAIN 

A2 

MAIN 

A3 

MAIN 

AA 

MAIN 

A5 

MAIN 

A6 

MAIN 

A7 

MAIN 

A8 

MAIN 

A9 

MAIN 

50 

MAIN 

51 

MAIN 

52 

MAIN 

53 

MAIN 

5A 

MAIN 

55 

MAIN 

56 

MAIN 

57 

MAIN 

58 

MAIN 

59 

MAIN 

60 

MAIN 

61 

MAIN 

62 

MAIN 

63 

MAIN 

6A 

I 


Li  t  ‘V _ 1 _ 1  _ _ _ _  ..  J J  Id _ 


J 


WRITE<6*5)  TITLE*  LL 

MAIN 

63 

WR1T6<6,15) 

MAIN 

66 

READ (5^ 20)  C*  D,  H*  RHOA*  RHOB#  RHOC#  RHOAP*  RHOBP#  RHUCP*  G 

MAIN 

67 

READ(5»20)  ( GG(1) #I>1» 2) 

MAIN 

69 

REA0(3*20)  ( RRHO ( I  )* 1*1*  9) 

MAIN 

70 

DO  600  IG  -  1*2 

MAIN 

71 

G  •  GG(IG) 

MAIN 

72 

DO  670  ID  ■  1*4 

MAIN 

73 

CALL  AXIS (0.0* 0.0* 4HC+CP»-4* 10.0*0.0*0. 0*0. 2) 

CORRB 

1 

CALL  AXIS(0.C*0.0»2HPP*2i 10.0*90.0*0.0*0.2) 

MAIN 

75 

GAP  »  FLOAT (ID) 

MAIN 

76 

DO  660  IRHO  -  1*9 

MAIN 

7B 

RHOC  -  RRHO(IRHO) 

MAIN 

79 

RHOCP  ■  RHOC 

CORRB 

2 

IPCC  •  0 

MAIN 

80 

DO  659  ICC  ■  1»41 

MAIN 

B1 

C  -  (ICC-1)  *  2.0  /  40.0 

MAIN 

82 

IF(C  .EQ.  0.0)  C  -  0.001 

MAIN 

83 

IF ( C  .GT.  GAP)  GO  TO  699 

MAIN 

84 

C  ■  C  /  2.0  t  CP  ■  C 

CORRB 

3 

D  ■  GAP  ••  C  -  CP 

CORRB 

4 

H  ■  0.0  %  HP  ■  0.0 

CORRB 

9 

RHOA  -  0.0  $  RHOAP  •  0.0  1  RHOB  •  0.0  1  RHOBP  -  0.0 

MAIN 

86 

IP ( G  .EQ.  0.209)  HF  -  9.762E9 

MAIN 

87 

I F ( G  .EQ.  0.290)  HF  -  3.064E9 

MAIN 

88 

IP  ( G  .EQ.  1.0)'  HF  >  3.333E9 

MAIN 

69 

IF ( G  .EQ.  0.260)  HF  -  3.33E9 

MAIN 

90 

IP(IC  .EQ.  2)  HP  ■  H 

MAIN 

91 

CAPL  *  RADIUS  +H+A+B4C+D+AP+BP+CP+HP 

MAIN 

93 

Cl  ■  (VT-VB)  /  (AL*VB) 

MAIN 

96 

P  ■  0.0 

MAIN 

97 

DX,  ■  AL  /  100.0 

MAIN 

98 

X  -  DX  /  2.0 

MAIN 

99 

DO  100  I  ■  1*99 

MAIN 

100 

IP ( RHOZ  .EQ.  0.0  .OR.  RHOT  . EQ.  0.0)  GO  TO  190 

MAIN 

101 

P  ■  P  +  ( SORT (RHOZ/RHOT) )  *  DX  *  i 1. 0+ (C1*CAPL+C 1 *P ) / < 1 .0+Cl*X ) ) ** 

MAIN 

102 

10.9/U.O+ALPHA+CAPL/E)  *♦  0.9 

MAIN 

103 

X  ■  X  +  DX 

MAIN 

104 

100  CONTINUE 

MAIN 

103 

c 

190  WRITE ( 6* 40)  P*  Cl*  BETA 

MAIN 

106 

130  CONTINUE 

MAIN 

107 

CAPL1  ■  H  +  A  +  C  /  2.0 

MAIN 

108 

PA  ■  0.0 

MAIN 

109 

DX  -  AL  /  100. 0 

MAIN 

110 

X  ■  AL  -  DX  t  2.0 

MAIN 

111 

DO  200  I  ■  1*99 

MAIN 

112 

IFtRHOZ  .IQ.  0.0  .OR.  RHOA  .EQ.  0.0)  GO  TO  210 

MAIN 

113 

PA  -  PA  +  ( SORT ( RHOZ /RHOA) ) 40X4(1 .0+ ( C14H+C14PA ) / (1 , 0+C IPX  )  )**0, 9/ 

MAIN 

114 

1  ( 1 .0  + AL PHA*H/E )  **  0.9 

MAIN 

113 

X  «  X  -  DX 

MAIN 

116 

IF(PA  .GE.  A)  GO  TO  210 

MAIN 

117 

200  CONTINUE 

MAIN 

1  IB 

210  DELLA  «  AL  -  X 

MAIN 

119 

PC  •  0.0 

MAIN 

120 

BETA  -  1.0 

MAIN 

121 

CAPVS  -  9.0E9 

MAIN 

122 

DX  •  (Au-DELLA)  /  100.0 

MAIN 

123 

X  -  (AL-DELLA)  -  DX  t  2.0 

MAIN 

124 

DO  220  1  ■  1*99 

MAIN 

129 

IF ( RHOC  .EQ.  0.0)  GO  TO  1090 

MAIN 

1 26 

c 

MAIN 

127 

c 

NEGATIVE  BETA  FOR  DETONATION  LAYER. 

MAIN 

128 

c 

♦44  INCLUDES  PRESSURE  DROP  WITH  EXPANSION.  ^ 

MAIN 

129 

c 

MAIN 

130 

IF  t  G  .IT.  0.0)  BETA  ■  1.0  -  G  *  (C /C C  +  PC* VDE T/C ARVS ) )  1  <0.S*RHQC* 
1RH0Z*VT*VT/!  r$QRT(RH0C>+SQRT(RH0l))**2n 
IF ( BET  A  . LT,  0.0)  BETA  •  0.0 
IF  C  G  .LT.  0.0)  60  TO  1050 

RHOJ  «  RHOZ  *  E  /  ((E+ALPHA*(H  +  A4C))*(1.04C1*(H+A*I*C)/(1.0+C1*X))) 


.LT.  0.0)  60 
■  RHOZ  *  E  / 
RHOJ  f  RHOC 
VB  '+  (VT-VB) 
5  ■  VJ  /  (1.0 
:K  ■  2.0  *  G 


-VB)  ♦  X  /  AL 
CAP  VS  ■  VJ  /  (1. O+SQRT (RHOC /RHOJ  ) ) 

DGREEK  ■  2.0  *  6  *  HF  /  CAPVS**2 
AKK  ■  1.0 

IF  (  G  .GT.  0.0)  AKK  -  1.0  f  ( 1 .0+ ( 1 ,0-G >*C Vf ( G*C VH20)  ) 

:  IF ( G  »EQ .  1.0)  60  TO  1000 

TEMP  ■  ( 0 . 5* RHOJ  *RH0C*V J*V  J/ ( 2. 0* (SORT (RHOJ ) +IQRT (RHOC  )  )**2 ) )  /  (< 
11.0-G) *CV) 

IF (TEMP  .LT.  373.0)  60  TO  1000 
GO  TO  1010 
1000  TEMP  ■  373.0 
1010  UORKEF  ■  <  TEMP-373,0)  /  TEMP 
GAM  «  WORKEF  *  AKK 

TF AC  -  1.0  -  (GAM/RO)  *  (l.O-GAM/RO)  A  (1.0-D6REEK) 

IF ( TFAC  .LE.  0.0)  GO  TO  1030 
GO  TO  ICAO 
1030  BETA  •  1.0 
GO  TO  1050 

1040  V  ■  (1. 0-3QR  T (TF  AC ) )  /  (l.O-GAM/RO) 

FRHO  -  SORT(RO) 

FRHO  •  (FRHO*FRHO+FRHO-1 .0 )  /  (l.O+FRHO) 

BETA  -  1.0  1  (1.  O+GAMM  V*V *FRHD-D GREEK ) ) 

IF ( BETA  .LT.  0.0)  BETA  -  0.0 
103u  RH01'  ■  BETA  *  RHOC 

IF ( RH01  .LT.  RHOC)  RH01  ■  RHOC 

IF ( RHOZ  .EQ.  0.0  .OR.  RH01  .EQ.  0.0)  GO  TO  230 

PC  •  PC  +  ( S QR T ( RH0Z/RH01)  )  *  OX  *  ( 1.0+( Cl* ( H+A+PC ) ) / (1. 0+Cl*X )) 
1**0. 3  f  (1.04ALPHA*(H4A)/E )  **  0.3 
X  ■  X  -  DX 

IF ( PC  .GE.  C)  GO  TO  230 
220  CONTINUE 

230  DELLC  «  AL  -  DELLA  -  X 
P9  ■  0.0 

DX  ■  ( AL-OELL  A-DELLC )  /  100.0 
X  -  (AL-DiLLA-DELLC)  -  OX  t  2.0 
DO  240  I  -  1,99 

IF ( RHOZ  .EQ.  0.0  .OR.  RH03  .EQ.  0,0)  GO  TO  250 

PB  ■  PB  +  ( S QRT ( RHOZ/RHOB )  )  *  OX  *  ( 1 . 0+ ( Cl* ( H  + A +C  +  P B ) ) / ( 1 .04C1* X ) 
1 )  **0«  5  t  1 1 , 0+AL  PH A* ( H+A+C ) /E )  **  0.3 
X  ■  X  -  DX 

IF ( PB  .GE.  B)  GO  TO  230 
240  CONTINUE 

230  DELLB  «  AL  -  DELLA  -  DELLC  -  X 
DELL  1  ■  DELLA  +  DELLB  +  OE  LLC 
C  WR'ITE  ( 6, 92  )  0ELL1 
C  WRITE!  6,  B  2  )  DELLA,  DELLB,  DELLC 

RH031  ■  RHOZ  *  E  I  ( ( 1 . 0+Cl*(  H+A+B+C  ) '/  ( 1. 0+Cl*(  AL-DE  LL1 ) ) )  *  (E  +  AL 
1PHA*( AL-DELL 1 ) ) ) 

El  "  E  f  (1.0+AL PHA*(H+A+9 +C ) /E) 

QCOL  ■  AK  *  El  *  RH0J1 
RCOL  •  RHOB  *  B 
0TS3  ■  B 
DIB  ■  B 
RHOP  -  RHOB 

IF ( QCOL  .LE.  RCOL)  GO  TO  260 

RCOL  ■  RCOL  +  RHOC  *  C 

DISS  ■  DISS  +  C  a c 

DIB  •  C 


(1.0+(1.0-G)*CVf(G*CVH20) > 


/  100,0 
■  OX  t  2.0 


MAIN 

131 

MAIN 

132 

MAIN 

133 

MAIN 

134 

MAIN 

135 

MAIN 

136 

MAIN 

137 

MAIN 

136 

MAIN 

139 

MAIN 

140 

MAIN 

141 

MAIN 

142 

MAIN 

143 

MAIN 

144 

MAIN 

145 

MAIN 

146 

MAIN 

147 

MAIN 

148 

MAIN 

149 

MAIN 

190 

MAIN 

131 

MAIN 

152 

MAIN 

193 

MAIN 

134 

MAIN 

155 

MAIN 

136 

MAIN 

137 

MAIN 

13B 

MAIN 

119 

MAIN 

160 

MAIN 

161 

MAIN 

162 

MAIN 

163 

MAIN 

164 

MAIN 

165 

MAIN 

166 

MAIN 

167 

MAIN 

16B 

MAIN 

169 

MAIN 

170 

MAIN 

171 

MAIN 

172 

MAIN 

173 

MAIN 

174 

MAIN 

175 

MAIN 

176 

MAIN 

177 

MAIN 

178 

MAIN 

179 

MAIN 

1B0 

MAIN 

1  81 

MAIN 

1B2 

MAIN 

163 

MAIN 

1*4 

MAIN 

1 85 

MAIN 

166 

MAIN 

187 

MAIN 

168 

MAIN 

189 

MAIN 

190 

MAIN 

191 

MAIN 

192 

MAIN 

193 

MAIN 

194 

45 


RHOP  ■  RHOC 

HAIN 

193 

IF ( OCOL  . LE.  RCOL  )  60  TO  260 

HATH 

196 

RCOL  ■  RCOL  t  RHOA  ♦  A 

HAIN 

197 

DOSS  -  OISS  +  A 

HAIN 

198 

DIS  •  A 

HAIN 

199 

RHOP  ■  RHOm 

HAIN 

200 

X  F ( QCOL  .10.  RCOL  )  GO  TO  260 

HAIN 

201 

X F ( A+B+C  .EQ.  0.0)  GO  TO  2  53 

HAON 

202 

RHOP  -  RCOL  /  (A+B+l) 

HAIN 

203 

GO  TO  255 

HATN 

204 

253 

RHOP  «  0.0 

HAON 

205 

255 

GO  TO  270 

HAIN 

2  06 

260 

IP ( RCOL  .EQ.  0.0)  GO  TO  263 

HAIN 

207 

FF  «  ( QCOL-RCOL )  t  RCOL 

HAIN 

206 

GO  TO  267 

HATN 

209 

265 

FF  -  0.0 

HAIN 

210 

267 

RCOL  ■  RCOL  -  RHOP  ♦  013  *  FF 

HAIN 

211 

DOSS  •  DISS  -  DIS  ♦  FF 

HAIN 

212 

IF ( DIS S  .EQ.  0.0)  GO  TO  266 

HAIN 

213 

RHOP  -  RCOL  f  OISS 

HAIN 

214 

GO  TO  270 

HAIN 

213 

268 

RHOP  -  0.0 

HAIN 

216 

270 

CONTINUE 

HAIN 

217 

VI  •  VB  ♦  (VT-VB)  ♦  (1.0-DELL1/AL) 

MAIN 

21B 

PRES  -  0.3  ♦  RHOP  *  RH0J1  *  VI  ♦  VI  /  < SQRT< RHOP ) 6SQRT ( RH0J1 )) t*2 

MAIN 

219 

c 

WRITE ( 6* 30 )  CAPLf  PA*  PB,  PC ,  Vli  RHOP,  RH01*  RH0J1*  PRES 

HAIN 

220 

F  •  1,0  f  (1.0+(SQRT(PRES7RHOJ1) ) M ALPHA* VI) ) 

MAIN 

221 

1F( (A6RHOA+B6RHOB+C4RHOU  .EQ.  0.0)  GO  TO  274 

HAIN 

222 

Y1  ■  AK  *  El  +  RH0U1  /  ( 1. 0+E 1/ ( A6RH0A+BPRH0B+CPRHQ1 ) ) 

MAIN 

223 

GO  TO  276 

MAIN 

224 

274 

Yl  -  0.0 

MAIN 

223 

276 

DELY1  ■  -Yl  ♦  (l,0-SQRT((1.0+ALPHA*D/E)/< 1 .0+ALPHA*F *D7£ ) )  ) 

MAIN 

226 

P2  -  0.0 

MAIN 

227 

CAPL2  -  H+  A  +  B  +  C+  D 

HAIN 

2  2B 

DX  ■  ( AL-DELL1)  /  100.0 

MAIN 

229 

X  ■  (AL-DELL1)  -  DX  /  2.0 

MAIN 

2  30 

DO  300  I  -  1*09 

MAIN 

231 

P2  -  P2  +  OX  ♦  <1.0+C16CAPL2/U.0+C1*X)) 

MAIN 

232 

X  -  X  -  DX 

MAIN 

233 

IF(P2  .GE.  DELY1 )  GO  TO  310 

MAIN 

234 

300 

CONTINUE 

MAIN 

2  35 

310 

DELL1P  ■  AL  -  DELL1  -  X 

MAIN 

2  36 

c 

WRITE ( 6* 60 )  F»  Yl,  DELY1,  P2,  DELL1P 

MAIN 

237 

c 

WR1TE( 6*62 )  RH02 

MAIN 

236 

PBP  «  0.0 

MAIN 

2  39 

OX  ■  ( AL-DELL1-DELL1P)  /  100.0 

MAIN 

240 

X  ■  (AL-DELL 1-DE  LL1P )  -  DX  t  2.0 

MAIN 

241 

DO  400  I  -  1,99 

HAIN 

242 

OF ( RHOZ  , EQ.  0.0  .OR.  RHOBP  .EQ.  0.0)  GO  TO  410 

MAIN 

243 

PBP  «  PBP  +  (  S Q R T  ( RHQZ / RHOBP )  )  ♦  OX  ♦  (l.  0+<ClMC  APL2+PBP)  >/(1.04C 

MAIN 

244 

n*xn  **  0.3  /  (i.o+alpha*capl2/ei  **  o.s 

HATN 

245 

X  -  X  -  DX 

MAIN 

246 

OF (PBP  .GE.  BP)  GO  TO  410 

MAIN 

2  47 

400 

CONTINUE 

MAIN 

24B 

410 

DELLBP  •  AL  -  DELL1  -  DELL  IP  -  X 

MAIN 

249 

c 

WROTE (6*70)  PBP*  BP,  DELLBP 

HAIN 

230 

PCP  ■  0.0 

MAIN 

251 

BETA  -  1,0 

HAIN 

2  32 

DX  -  (  AL-DELL 1-DELL  IP-DELL  BP )  t  100.0 

MAIN 

233 

X  ■  (AL-DELL l-DELL 1P-DELLBP)  -  DX  /  2.0 

HAIN 

234 

DO  420  I  ■  1*»  99 

MAIN 

255 

c 

MAIN 

2  56 

c 

♦  ♦♦ 

NEGATIVE  BETA  FOR  DETONATION  LAVER. 

INCLUDES  PRESSURE  DROP  WITH  EXPANSION 

MAIN 

257 

c 

♦  64 

MAIN 

256 

c 


c 

MAIN 

2  59 

IF (RHOCP  .EQ.  0.0)  60  TO  2050 

MAIN 

2  60 

IF  ( G  .LT.  0.01  B'ETA  -  1.0  -  G  4  (C  P/t  C P+PCP4VDET/CAPVS ) )  t  (0.S4 

MAIN 

261 

lRH0C?4RH0Z4WT4VT/(  < SQRTtRHOC P ) +SQRT( RHOZ > )442 ) ) 

MAIN 

2  62 

IF(BETA  .LT.  0.0)  BETA  •  0.0 

MAIN 

263 

IF ( G  .LT.  0.0)  60  TO  2030 

MAIN 

264 

RHOd  ■  RH0Z  ♦  E  /  <<E+ALPHA4(CAPL24BP> >4(1.0+C14(CAPL2+BP4PCP>/(1. 

MAIN 

2  63 

104C14X) ) ) 

MAIN 

2  66 

RO  •  RHOJ  /  RHOCP 

MAIN 

267 

VU  •  VB  +  (VT-VB)  4  X  /  AL 

MAIN 

268 

CAPVS  ■  VJ  /  (1.0+SQRT(RH0CP/RHCU)> 

MAIN 

269 

DGREEK  -  2.0  4  G  4  HF  /  C APVS**2 

MAIN 

2  70 

AKK  ■  1.0 

MAIN 

271 

IF (S  . 6T .  0.0)  AKK  -  1.0  /  ( 1. 0  +  ( 1 .0-6 )4C V/ ( G4C VH20) ) 

MAIN 

272 

C 

I F  ( G  ,EQ.  1.0)  60  TO  2000 

MAIN 

2  73 

TEMP  ■  (0.34RHO|4RHOCP4VJ4VJ/(2.04<SORT(RHQJ)+SCIRT(RHOCP))442))  / 

MAIN 

2  74 

i  (d.o-e)*cv> 

MAIN 

2  73 

T F  t TEMP  .LT.  373.0)  60  TO  2000 

MAIN 

276 

GO  TO  2010 

MAIN 

2  77 

2000 

TEMP  ■  373.0 

MAIN 

2  7B 

2010 

WORKEF  -  (TEMP-373.0)  /  TEMP 

MAIN 

2  79 

GAM  •  WORKEF  4  AKK 

MAIN 

2  80 

TF AC  -  1.0  -  (GAM/RO)  4  U.O-GAM/RO)  4  (1.0-DGREEK) 

MAIN 

2  81 

IF  ( TF AC  .LE.  0.0)  60  TO  2030 

MAIN 

2  62 

GO  TO  2040 

MAIN 

2  83 

2030 

BETA  -  1.0 

MAIN 

2  84 

GO  TO  2030 

MAIN 

2  B3 

rv 

o 

-r* 

o 

V  ■  (l.O-SQRT(TPAC))  /  (l.O-GAM/RO) 

MAIN 

286 

FRHO  -  BQRT(RO) 

MAIN 

287 

FRHO  -  (FRH04FRHQ+FRH0-1 .0)  /  (l.O+FRHO) 

MAIN 

288 

BETA  ■  1.0  /  U.O+GAM4(V4V4FRHQ-OGREEK)) 

MAIN 

2  89 

IF ( BETA  .LT.  0.0)  BETA  ■  0.0 

MAIN 

2  90 

2050 

RH02  •  BETA  4  RHOCP 

MAIN 

2  91 

IF ( RH02  .LT.  RHOCP)  RHQ2  ■  RHOCP 

MAIN 

2  92 

I F ( RHOZ  .EQ.  0.0  .OR.  RH02  .EQ.  0.0)  GO  TO  430 

MAIN 

2  93 

PCP  »  PCP  +  ( SORT ( RHOZ/RH02 )  )  *  OX  4  (1.0+C14(CAPL2+BP+PCP )  /  (1.0 

MAIN 

2  94 

1+C14X))  44  0.3  /  ( 1.0+ALPHA4(CAPL2+BF)/E)  44  0.5 

MAIN 

2  93 

X  ■  X  -  DX 

MAIN 

296 

IF ( PC P  .GE.  CP)  GO  TO  430 

MAIN 

2  97 

4  20 

CONTINUE 

MAIN 

2  96 

430 

OELLCP  ■  AL  -  DELL1  -  0ELL1P  -  DELLBP  -  X 

MAIN 

2  99 

C 

WRITE  (6. 94)  PCP..  CP.  OELLCP 

MAIN 

3  00 

PAP  ■  0.0 

MAIN 

301 

DX  •  l AL-DELL1-DELL1P-DELLBP-DELLC  P)  /  100.0 

MAIN 

302 

X  -  ( AL-OELL 1-DELL 1P-0ELLB P-OE LLC P )  -  OX  /  2.0 

MAIN 

303 

DO  440  I  ■  1.99 

MAIN 

3  04 

IF ( RHOZ  .EQ.  0.0  .OR.  RHOAP  .EQ.  0.0)  GO  TO  430 

MAIN 

303 

PAP  -  PAP  +  ( SORT ( RHOZ /RHOAP ) )  4  OX  4  ( 1 . 0+C 14 ( C APL2+BP+CP +PAP ) / ( 1 

MAIN 

306 

1.0+C14X))  44  0.3  /  (1.0+ALPHA4(CAPL2+BP+CP)/E)  44  0.3 

MAIN 

3  07 

■'  X  «  X  -  DX 

MAIN 

3  08 

I F ( P AP  • GE .  AP)  GO  TO  430 

MAIN 

3  09 

440 

CONTINUE 

MAIN 

310 

450 

OELLAP  ■  AL  -  DELL1  -  0ELL1P  -  DELLBP  -  OELLCP  -  X 

MAIN 

311 

DELL 2  «  DELLAP  +  DELLBP  ♦  DELLCP 

MAIN 

312 

C 

WRITE (6.96)  PAP.  AP>  DELL APy  DELL2 

MAIN 

313 

V2  •  VB  +  (VT-VB)  4  ( AL-DE LL1-0ELL 1P-DELL2 )  /  AL 

MAIN 

314 

IF ( AP  +  BP  +  C P  .EQ.  0.0)  GO  TO  431 

MAIN 

315 

RHOPP  •  ( AP4RH0AP+BP4RH0BP+CP4RHUCP)  /  (AP+BP+CP) 

MAIN 

316 

GO  TO  432 

MAIN 

317 

451 

RHOPP  «  0.0 

MAIM 

316 

432 

RH0J2  «  RHOZ  4  E  /  ( ( ( 1.0+C14( CAP12  +  AP+BP+CP ) /(l. 0+C 14( Al-DELll- 

MAIN 

319 

1DELL1P-DELL2) ) )4(E+ALPHA4(  Al-DELL 1-DELL1 P-DE  LL2 ) ) ) ) 

MAIN 

3  20 

E 2  «  E  '/  (1.0  +  ALPHA4(CAPL2+AP  +  BP+CP)/E‘  .7 

QCOL  •  AK  4  >E2  4  RH0J2 

MAIN 

321 

MAIN 

3  22 

RCOL  -  RHOBP  4  AP 

MAIN 

323 

OISS  *  AP 

MAIN 

324 

DIS  ■  AP 

MAIN 

325 

RHOPP  -  RHOAP 

MAIN 

326 

IP ( QCOL  . L E.  RCOL )  GO  TO  460 

1 

MAIN 

327 

RCOL  ■  RCOL  4  RHQCP  *  CP 

MAIN 

328 

OISS  •  OISS  +  CP 

MAIN 

329 

OIS  ■  CP 

MAIN 

330 

RHOPP  •  RHOCP 

MAIN 

331 

IF ( QCOL  .LE.  RCOL)  GO  TO  460 

MAIN 

332 

RCOL  •  RCOL  4  RHOAP  *  AP 

MAIN 

333 

OTSS  •  OISS  +  AP 

MAIN 

3  34 

DIS  •  AP 

MAIN 

333 

RHOPP  •  RHOBP 

MAIN 

336 

IF(QCOL  .LE.  RCOL)  GO  TO  460 

MAIN 

3  37 

IP ( AP+BP+CP  .EQ.  0.0)  GO  TO  4B3 

MAIN 

338 

RHOPP  ■  RCOL  /  ( AP+BP+CP) 

MAIN 

339 

GO  TO  4*3 

MAIN 

340 

453 

RHOP  -  0.0 

MAIN 

341 

453 

GO  TD  470 

MAIN 

342 

460 

IP ( RCOL  .EQ.  0.0)  GO  TO  463 

MAIN 

343 

FP  •  (QCOL -RCOL)  t  RCOL 

MAIN 

344 

GO  TO  467 

MAIN 

343 

463 

FF  •  0.0 

MAIN 

346 

467 

RCOL  •  RCOL  “  RHOPP  4  FF  4  OIS 

MAIN 

347 

DISS  ■  DISS  -  OIS  4  FP 

MAIN 

348 

IFIDISS  .EQ.  0.0)  GO  TO  46« 

MAIN 

349 

RHOPP  ■  RCOL  /  OISS 

MAIN 

310 

GO  TO  470 

MAIN 

331 

468 

RHOP  ■  0.0 

MAIN 

352 

470 

CONTINUE 

MAIN 

333 

PRES  ■  0.3  4RHOPP4RHO)IZ4V2*V2/(SORT(  RHOPP)  4SQRT( RHOJ  2  ) )  442 

MAIN 

394 

WRITE ( 6*  98 )  RHOPP.  RH0J2*  V2 

MAIN 

333 

FP  •  1,0  t  ( 1.04 (SORT( PRE5/RH0J2) )/(AL PHA*V2) ) 

MAIN 

3  36 

IF( ( AP4RHOAP+BP4RHOBP4CP4RHOCP)  .EQ.  0. 

0)  GO  TO 

474 

MAIN 

337 

Y2  ■  AK  *  E2  /  (1.04E24RHOJ2/UP4RHOAP4BP4RHOBP4CP4RKOCP)) 

MAIN 

338 

GO  TO  476 

MAIN 

339 

474 

Y2  ■  0.0 

MAIN 

360 

476 

DEL Y2  -  -Y2  4  (1,0-SQRT((1.0+AL PHA4HP/E) / ( 1 , 04ALPHA4FP4HP/E ) ) ) 

MAIN 

361 

P4  •  0.0 

MAIN 

362 

DX  ■  UL-(DELL140ELL1P4DEL12))  //  100.0 

MAIN 

363 

X  ■  99.3  4  DX 

MAIN 

364 

DO  300  I  -  1 » 99 

MAIN 

363 

P4  ■  P4  4  DX  4  (1.04C14CAPL/U.04C14X)) 

MAIN 

3  66 

X  »  X  -  DX 

MAIN 

367 

IF ( P4  .GE.  DEL Y2 )  GO  TO  310 

MAIN 

3  68 

300 

CONTINUE 

MAIN 

3  69 

310 

DELL2P  ■  AL  -  (DELL14DELL1P4DELL2)  -  X 

MAIN 

3  70 

WRITE ( 6v  80 )  Y2»  V2»  PRES*  F t  DELY2*  P4* 

DELL2P* 

FP 

MAIN 

371 

PP  •  0,0 

MAIN 

3  72 

DX  ■  AL  -  (DELL14DELL&P4DELL24DELL2P) 

MAIN 

373 

OX  ■  DX  /  100,0 

MAIN 

374 

X  -  DX  /  2.0 

MAIN 

373 

DO  600  I  -  1*99 

MAIN 

3  76 

TF<  RHOZ  .EQ.  0.0  .OR.  RHOT  .EQ.  0.0)  GO 

TO  630 

MAIN 

3  77 

PP  ■  PP  4  (SORT ( RHQZ/RHOT) )  4  OX  4  1 1,04 ( C14C APL4C14PP )/( 1 .04C14X) 

MAIN 

378 

1) 440. 3/(1.04ALPHA4CAPL/E>  **  0.3 

MAIN 

379 

X  -  X  4  DX 

MAIN 

380 

600 

CONTINUE 

MAIN 

381 

PP IN  •  PP  /  2.34 

MAIN 

382 

630 

CONTINUE 

CORRB 

6 

IF(C  .GT.  2.0  .OR,  PP  .GT.  2.0)  GO  TO  655 

MAIN 

384 

IPCC  ■  IPCC  4  1 

PLC(ICC)  ■  C  4  CP  %  PLPPUCC)  «  PP 

MAIN 

383 

CORRB 

7 

659 

CONTINUE 

MAIN 

380 

ncapccui  ■  o.o 

MAIN 

389 

puPMiPCc+n  «  o.o 

MAIN 

3  90 

PLCUPCCA2)  •  0.2 

MAIN 

391 

PLPTMIPCC  +  2)  -  0.2 

MAIN 

392 

CALL  HNE(PLt*PLP°»IPCCil»0>0) 

MAIN 

393 

660 

CONTINUE 

MAIN 

394 

CALL  PLTP6E 

MAIN 

395 

670 

CONTINUE 

MAIN 

396 

680 

CONTINUE 

MAIN 

397 

700 

CONTINUE 

MAIN 

390 

3000 

CONTINUE 

MAIN 

399 

STOP 

MAIN 

4  00 

END 

MAIN 

4  01 
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